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Abstract 

The “simplified shear solution” method is presented for approximating the through-thickness shear 
stress distribution within a composite laminate based on laminated beam theory. The method does not 
consider the solution of a particular boundary value problem, rather it requires only knowledge of the 
global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to 
lamination theory in that ply level stresses can be efficiently determined from global load resultants (as 
determined, for instance, by finite element analysis) at a given location in a structure and used to evaluate 
the margin of safety on a ply by ply basis. The simplified shear solution stress distribution is zero at free 
surfaces, continuous at ply boundaries, and integrates to the applied shear load. Comparisons to existing 
theories are made for a variety of laminates, and design examples are provided illustrating the use of the 
method for determining through-thickness shear stress margins in several types of composite panels and 
in the context of a finite element structural analysis. 


Introduction 

It is known that by employing lamination theory, one can determine the in-plane stress distribution in 
each layer of a laminate from the knowledge of the applied force and moment resultants. This information 
enables determination of ply-by-ply margins for the laminate, which are needed for design and sizing. 

The distribution of the interlaminar shear stresses in each layer from the applied shear resultants are not 
readily available from the standard lamination theory equations. These stresses are also needed to enable 
design and sizing of laminates subjected to global shear loads. The present report provides an analytical 
method, based on laminated beam theory with shear loading, for determining the interlaminar shear stress 
distribution in a laminate from a given applied shear resultant. It should be emphasized that the classical 
or the higher-order (e.g., first-order shear deformable) plate theory cannot be employed to determine the 
correct interlaminar shear stress distribution through the laminate thickness from the knowledge of the 
force, moment, and shear resultants. Indeed, the classical plate theory provides identically zero 
interlaminar shear stresses, whereas the higher-order plate theories provide piece-wise profiles that are 
discontinuous at the ply interfaces. It should also be noted that the simple shear solution provided herein 
does not involve the solution of a particular boundary value problem. Rather, it is only assumed that the 
force, moment, and shear resultants are known at a particular location in the laminate, and the solution is 
independent of the source of these known resultant quantities. Conversely, in order to obtain the 
interlaminar shear distribution via integration of the differential equations of equilibrium, the solution of a 
boundary value problem would be required to determine the variation of the stress fields in each ply. 
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The presented theory has been validated for laminates with one to five layers composed of both 
isotropic and anisotropic materials by comparison to approximate methods (HOTFMG (Aboudi et al., 
1999), HyperSizer Joints (Zhang et al., 2006)) and an exact plate solution (Williams, 1999). Finally, the 
theory is implemented for the practical analysis of composite stiffened panels in the FlyperSizer software 
(Collier Research Corp., 2007). Results, including predicted margins of safety, are given for a 
graphite/epoxy laminate, a honeycomb sandwich panel, an I-stiffened panel, and laminate components of 
a pressure shell structure. 


Formulation 


We seek to determine the distribution of the shear stress in a laminate given the global applied shear 
resultant, Q. Consider a laminate (fig. 1) with an applied shear force P consisting of N layers, with the 
layer number denoted by k. Note that the positive directionality of P is motivated by that which causes a 
positive internal shear stress (see app. B). Following the concept of a composite beam (Beer and 
Johnston, 1981), the actual laminate cross-section is replaced with an effective cross-section. As shown in 
figure 2, the original width of each layer ( b ) is transformed to b n k , where, 


E x is the Young’s modulus of layer k in the x-direction, and E x is the average Young’s modulus of the 
laminate in the x-direction. E x serves only as a normalization factor which allows the ply stiffnesses to 
be scaled with respect to each other. Alternatively, the Y oung's modulus of one particular ply may be 
used. Note that in bending, under a plane strain condition, where, £ v =0, E x should be modified to 
E x —> E x /[l - v X yVy X ] (see app. A). The neutral axis for the laminate is located at the centroid of the 
transformed laminate cross-section and is given (with respect to the laminate midplane coordinate system) 

by, 



( 2 ) 


where A k and 4 are the area and thickness of layer k, respectively. The distance of a point from the neutral 
axis is denoted by z = z-z* . 


Z 

vertical cut 

G x 



Figure 1. — Portion of a laminate subjected to a shear force, P, and the 
resulting normal stress distribution on an arbitrary cross-section. 
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Figure 2. — Original (actual) and transformed laminate cross-sections. 

The applied shear loading causes a normal stress distribution at an arbitrary vertical cross-section 
given according to beam theory (Beer and Johnston, 1981) by, 

Mz 



(3) 


where / is the moment of inertia of the transformed laminate cross-section (fig. 2) and M is moment at the 
vertical cross-section. Because a counter clockwise moment is considered positive, equilibrium dictates 
that M = -Px, thus, 


= 


Pxz 

I 


(4) 


To obtain the shear force, denoted by H, acting on an arbitrary horizontal cross-section at location 
z = z c , within ply number k (see fig. 3), consider the equilibrium of forces in the x-direction, which 
requires, 


-H + \a x dA = 0 (5) 

where dA is the area of the differential element indicated in figure 3. Note that the positive H sign 
convention indicated in figure 3 is consistent with that employed by Beer and Johnston (1981), see 
appendix B. Substituting eq. (4) into eq. (5) yields, 



( 6 ) 


where c - z N denotes the distance from the neutral axis to the bottom surface of the laminate (see fig. 3). 
Using dA - bn^dz , we obtain, 



n k z dz 


The integral appearing in eq. (7), which is denoted by Q'{k, z ) , can be decomposed as follows, 


( 7 ) 
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( 8 ) 


*N - 1 

Q'(k,z) = | bn k zdz = ^bn k zdz + J bn k+1 zdz + ...+ J bn N zdz 


z k + 1 


It follows that 


Q'(k,z) 



(9) 


Zc 


Zk 

/\ 

Z N = -C 


H 


neutral axis 


b n. 


Layer k 




— \ 



horizontal cut 


=r 


dz 

dA 


cross-section 

vertical cut 

Figure. 3. — Section of laminate acted upon by the horizontal shear force, H, and the distributed normal stress, a x . 
where z N =-c. Hence, the shear force, H, is given by, 


H = —Q'{k,z) 

The shear stress, x, is given by, 

__H 

xb 


( 10 ) 


( 11 ) 


where the negative sign is due to the opposite directionality of the internal shear force compared 

to the shear stress (see app. B). The applied shear resultant, Q, is related to the applied shear force, P, by 

Q = P/b , and it follows that, 


x (k,z) = 



( 12 ) 


In employing this equation, the moment of inertia of the transformed beam cross-section is given by, 


'=1 


—bn k t k +bn k t k 


f Z k +Z k-l^ 

2 

{ 2 J 



(13) 
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Thus, given eqs. (9) and (13), the shear stress distribution throughout the laminate can be determined 
from eq. (12). 

It should be emphasized that at the top and bottom of the laminate the shear stress must be equal to 
zero. In addition, the integral of t over the entire thickness of the laminate must be equal to the applied 
shear resultant, Q, that is, 


Z N 

| t dz = Q 

Zo 


(14) 


Verification Cases 

In order to verify the analytical formulation, we compare to results from the Higher-Order Theory for 
Functionally Graded Materials (HOTFGM) (Aboudi et al., 1999). The 2D HOTFGM problem domain 
considered consists of 20 x 200 subcells in the z- and x-directions, respectively, subjected to shear loading 
as shown in figure 4. It should be noted that the shear stresses shown in the results below were taken at 
x = 5 (middle of the domain) in order to ensure that the edge effects do not influence the results. 
Furthermore, comparison has been made to the HyperSizer Joints (Zhang et al., 2006) software, which is 
based on a cylindrical bending plate theory for adhesively bonded laminates. 

HOTFGM is a plane strain program, such that s y = 0. In the results cases 1 to 5, the Poisson ratio’s of 
the layers were assumed to be (nearly) identical, thus the plane strain condition coincides with the plane 
stress condition. In case 6, the effect of different layer Poisson ratios will be shown. 

Case 1: Monolithic Material 

In the case of a monolithic material, the shear stress distribution reduces to the beam theory parabolic 
expression. In this case, N= 1, the neutral axis is at the midplane (i.e., z* = 0 ), z 0 = c , and z 1 =—c, see 
figure 5. Therefore, from eq. (9), 



Since the cross-sectional area A = 2 be, /= b (2c) 3 /12, and Q = P/b, it follows that 


T 


3 P_ 
2 A 


f 

V 



(15) 


(16) 


which is the well-known shear stress distribution in a monolithic rectangular beam. 

The simplified shear solution (as implemented in FORTRAN) is compared with the HOTFGM 
solution and the HyperSizer Joints solution for the shear distribution in the monolithic material in 
figure 6. The three methods are in excellent agreement. 
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Figure 4. — HOTFGM representation of the laminate considered. 
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z 0 = c 
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z 1 = -C 

Figure 5. — Monolithic (single layer) laminate. 



Figure 6. — Monolithic material shear stress distribution (Case 1). 
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Figure 7. — Bi-material laminate shear stress distribution (Case 2). 


Case 2: Bi-Material (Stiffer on Top) 

The results for this case are shown in figure 7. The resulting shear distribution is piece-wise parabolic 
with the peak occurring in the stiffer material at the location of the neutral axis (z = 0.15). The results of 
the simplified shear solution, HOTFGM, and HyperSizer Joints coincide. 


Case 3: Bi-Material (Stiffer on Bottom) 

The results for this case are shown in figure 8. As in Case 2, the resulting shear distribution is 
piece-wise parabolic with the peak occurring in the stiffer material at the location of the neutral axis 
(z = -0.15). The results of the simplified shear solution, HOTFGM, and HyperSizer Joints are in 
excellent agreement. 


Case 4: Three Layers 

The results for this case are shown in figure 9. Because the laminate is symmetric, the shear 
distribution is symmetric as well, with maximum at the laminate midplane. The results of the simplified 
shear solution, HOTFGM, and HyperSizer Joints are in excellent agreement. 

Case 5: Four Layers 

The results for the four layer case are shown in figure 1 0. The peak occurs in the lower half of the 
laminate, which contains the stiffer materials. The results of the simplified shear solution, HOTFGM, and 
HyperSizer Joints are in excellent agreement. 

Case 6: Two Layers - Different Poisson Ratio’s 
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This case is identical to Case 3 except that the Poisson ratios of the layers have been changed. The 
bottom (stiffer) layer is given a Poisson ratio = 0.45, while the top (more compliant) layer is given a 
Poisson ratio = 0.05. The results shown in figure 1 1 indicate that the Poisson ratio mismatch causes a 
slight shifting of the shear stress distribution and amplification of the maximum compared to the Case 3 
results, which employed equal Poisson ratios. The simplified shear solution, HOTFGM, and HyperSizer 
Joints coincide. 



Figure 8. — Bi-material laminate shear stress distribution (Case 3). 
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In conclusion, it has been demonstrated that the present simplified shear solution coincides with the 
HOTFGM and HyperSizer Joints predictions for all of the above cases considered. 
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Cases 7 to 9: Comparison with Plate Theory 

Here we compare to results from Williams' (1999) multilength scale plate theory, which has been 
shown to coincide with Pagano's (1969) exact solution for plates in cylindrical bending, as well as 
HyperSizer Joints. It should be noted that Williams' plate theory considers sinusoidal pressure loading on 
the top of the plate as opposed to an applied shear load. The x x: distribution from a cross-section within 
the plate theory results were then integrated obtain the applied Q x used in the current simplified shear 
solution. This therefore mimics the HyperSizer situation where a Q x value is given for a point in the 
laminate with no knowledge of the global problem or boundary conditions used to generate the Q x load. 
In all three cases below, Q x = 0.32 N/mm. We consider three asymmetric laminates, a cross-ply, a 
[0°/9074579070°], and a [-3 0790745' 760730°], where the asymmetry arises due to different ply 
thicknesses, which are [0.2, 0.4, 0.2, 0.1, 0.1] mm for all three laminates. The material properties of the 
plies are given in table 1, where the standard expression for v vx is given in eq. (30) in appendix A. The 
resulting shear stress distributions for these laminates are shown in figures 12 to 14, respectively. 


TABLE 1— MATERIAL PROPERTIES FOR SIMPLIFIED SHEAR SOLUTION COMPARISON WITH PLATE THEORY. 


Ply angle, degree 

-30 

0 

30 

45 

60 

90 

E„ GPa 

2.192 

25 

2.192 

1.325 

1.068 

1 

E„ GPa 

1.068 

1 

1.068 

1.325 

2.192 

25 

Vxv 

0.4082 

0.25 

0.4082 

0.3245 

0.1989 

0.01 


Figure 12 shows that the shear stresses predicted for the cross-ply laminate coincide with the 
HyperSizer Joints solution, but deviate slightly from Williams' (1999) solution. Figure 13 shows that 
introducing a 45° ply results in some deviation in the top 90° and 45° plies. Introducing angle plies to the 
laminate configuration results in the distribution shown in figure 14. In conclusion, the simplified shear 
solution appears to give a good estimate of the shear distributions for the laminates considered, given the 
method's simplicity and efficiency. 



Figure 12. — Shear stress distribution for [0°/90 o /0 o /90 o /0°] laminate 
with ply thickness of [0.2, 0.4, 0.2, 0.1, 0.1] mm (Case 7). 
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Figure 13. — Shear stress distribution for [079074579070°] laminate 
with ply thickness of [0.2, 0.4, 0.2, 0.1, 0.1] mm (Case 8). 



Figure 14. — Shear stress distribution for [-30790745760730°] laminate 
with ply thickness of [0.2, 0.4, 0.2, 0.1, 0.1] mm (Case 9). 


Design Examples 

In order to design a composite laminate or stiffened panel that is subjected to external loading, the 
first step is to determine the induced local stresses in the structure. Then these stresses are compared to 
allowable stresses in order to determine the margin of safety of the considered composite panel. The 
propose simplified shear solution returns ply-level shear stresses at any point in the panel in global x, y, z 
laminate coordinates (see fig. 4). However, most failure analysis methods require stress allowables in 
terms of ply local (x\, x 2 , x 2 ) coordinates, where x 2 =z and x\ corresponds to the ply fiber direction. 
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Therefore, stresses in global coordinates obtained from the present analysis for any ply must be 
transformed into local ply coordinates according to: 

[113] cos0 sin0 

[T23 J -sinO cos0 

where 0 is the ply angle (between the fiber direction and the global x-direction). Note that in the previous 
sections, the shear stress t denotes either x x: , or x yz , depending on whether the applied shear force resultant 
Q stands for Q x or Q v (respectively). 

Once the ply local stresses are known, it is assumed that failure will occur when the following 
interaction equation is satisfied: 



( \2 

*13 

yFSU n . 


+ 


( _ 

*23 

Fsuji 


= 1 


08) 


where Fsun and Fsii 2 3 are the ply material shear stress allowables, respectively, which can be found in 
references such as M1L-HDBK-17. This interaction equation can be written as a margin of safety (see 
Zhang and Collier, 2004) as, 


MS = 


-1 


*13 


A2 


Fsu 


13 J 


*23 


A2 


Fsu 


23 J 


(19) 


where the margin of safety expresses how close to failure is the given panel. Thus, a positive margin of 
safety indicates that there is no failure for the considered loading, while a negative margin of safety 
indicates that failure has occurred. 

Finally, it should be mentioned that in the presence of combined applied shear loading (Q x and Q y ), 
the current analysis is applicable by superimposing the resulting stress distribution induced by each of Q x 
and Q y independently. 


Design Example 1 — Unstiffened Laminate 

This example considers an unstiffened graphite/epoxy laminate subjected to biaxial transverse shear 
loading (Q x = Q v = 200 lb/in.). This type of loading could result, for example, at the corner of a 
rectangular panel subjected to uniform pressure loading. A [45 7-45 7079070790745 7-45 °/0°] s layup is 
considered with uniform ply thicknesses of 0.005 in., thus forming a total thickness of 0.09 in. Ply 
material properties are given in table 2. 


TABLE 2.— GRAPHITE/EPOXY PLY MATERIAL PROPERTIES. 


Ex, 

e 2 , 

V 12 

FsUi 3 , 

Fsu 2 3 , 

Msi 

Msi 


ksi 

ksi 

23.35 

1.65 

0.32 

14.8 

5.32 


Figure 15 shows the resulting ply-by-ply shear stress distributions in the panel in the global 
coordinates as obtained by the simplified shear solution. Figure 16 shows the ply-by-ply stresses 
transformed from global coordinates to ply local coordinates. These shear stresses are then used within 
eq. (19) to calculate the margin of safety within each ply. Figure 17 shows these ply -by -ply margins of 
safety. Note that the minimum margin of safety is 0.03 which occurs in both plies 8 and 11. Since the 
margin of safety is positive but rather small, it can be concluded that the considered laminate is not highly 
over designed for the imposed loading. 
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Figure 15. — Shear stress distribution referred to the ply global coordinate system. 



Shear Stress (psi) 

Figure 16. — Shear stress distribution referred to the laminate local coordinate system. 
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Figure 17. — Margin of safety distribution within the considered laminate. 



Design Example 2 — Honeycomb Sandwich Panel 

This example considers a honeycomb sandwich panel, whose stress distribution is similar to an 
unstiffened composite laminate, except that the gradient of shear stress takes place almost entirely in the 
facesheets and is essentially constant across the core. The honeycomb sandwich facesheets considered are 
identical to the laminate considered in the previous example, with ply material properties given in table 2. 
The honeycomb core material properties are given in table 3, where it should be noted that the in-plane 
Y oung's modulus is taken to be extremely low. The shear strengths of the core differ in the x and y 
directions and are typical of Hexcel HRH 10 (Nomex/Phenolic) Honeycomb. The thickness of the core is 
1.5 in., and the applied loading is given by Q x = Q y = 1000 lb/in. 

The resulting stress distributions are shown in figure 1 8. The total thickness of the sandwich panel is 
1 .68 in., while the figure shows the distribution in the region 0.7 in. < z < 0.84 in., where z = 0 at the center of 
the sandwich. The shear stresses rise rapidly in the facesheet and are virtually constant in the core region. 
Referring to the local coordinates, the corresponding shear stress distributions are shown in figure 1 9. Once 
again, these local shear stresses are used to determine the margin of safety within the honeycomb sandwich, 
which is plotted in figure 20. As shown, the margins of safety in the relatively weak core (Fsuu = 0.11 ksi) 
are negative while margins of safety in the facesheets remain greater than 5. The negative margin of safety in 
the core indicates that core failure will occur under the applied loading, indicating the need for a stronger or 
thicker core. Conversely, the high margin of safety in the facesheet region indicates that the facesheet 
laminate is over designed for the applied shear loading. Note that sandwich panel facesheets are typically 
designed to carry in-plane loading, which is absent in the present example. 


TABLE 3 — HONEYCOMB CORE MATERIAL PROPERTIES 


T 1 

II 

Fsuxz, 

FSUy Z , 

psi 

ksi 

ksi 

0.0001 

0.225 

0.11 
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Figure 18. — Shear stress distribution referred to the sandwich panel global coordinate system. 



Shear Stress (psi) 

Figure 19. — Shear stress distribution referred to the sandwich panel local coordinate system. 
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Design Example 3 — I-Stiffened Panel 

This example investigates the shear stress distributions in an I-stiffened panel, as shown in figure 2 1 . 
If the shear force is aligned in the direction of the stiffeners (applied Q x ), then the stiffener webs 
contributed significantly to the transverse shear stiffness of the panel. On the other hand, if the shear force 
is applied transversely to the stiffeners (applied Q y ), the assumption is that the stiffener webs do not 
contribute to the shear stiffness of the panel. 

For the facesheets and flanges of the I-stiffeners, the laminate is treated in the same way as for 
unstiffened or sandwich panels by determining the Young's modulus and Poisson’s Ratio of each ply in 
the analysis (x or v) direction. For stiffener webs, if aligned with the analysis direction (for example, 
applied Q x on an I-Stiffened panel), the effective stiffness contribution of the web is determined by 
homogenizing the stiffener webs to an equivalent layer stiffness. This homogenization can be 
approximated by the simple expression: 


E eq = 


A“ F t 

yi \ 1 , web x , web 1 web 


( 20 ) 
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where A‘\ web is the “uncoupled” A u stiffness term of the web, E xweb is the effective elastic modulus of 

the web material, t weh is the thickness of the web, and S x is the stiffener spacing. A more extensive 
discussion of homogenization of the stiffener web can be found in Bednarcyk and Yarrington (2008). 

For shear force applied in the transverse direction ( Q v ), the stiffener web is assumed to be non- 
existent and the facesheet is treated as an unstiffened laminate. In the case of transverse shear force 
applied to a stiffened sandwich (such as a blade sandwich or an 1 - Sandwich, which incorporate bottom 
facesheets), the stiffener web stiffness is set to an extremely low value (to avoid numerical problems) and 
the panel is treated similarly to a honeycomb stiffened sandwich. It should be noted that the present 
simplified shear solution method provides a panel acreage analysis tool that captures interlaminar stresses 
away from stress concentrations, such as those that occur at free edges. 

This example considers a facesheet with ply material properties as given in table 2. The facesheet 
layup is, as previously considered, [45/-45/0/90/0/90/45/-45/0] s , with ply thicknesses of 0.005 in. thus 
forming a total thickness of 0.09 in. The upper flange and web of the I-stiffeners consist of [0/90/45/-45] s 
laminates with a ply thic kn ess of 0.00375 in., thus forming a total thic kn ess of 0.03 in. As for the lower 
flange, the layup is [0/90/45/-45/0/0] s , with ply thickness of [0.005/0. 005/0. 005/0. 005/0. 00375/0. 00375] s 
in., thus forming a total thic kn ess of 0.05 in. The applied shear loading is given by Q x = Q y = 200 lb/in. 

Stress profiles at cut A-A are shown in figures 23 to 25. Figure 23 shows the x xz and x yz shear stresses 
in global coordinates through the entire panel depth. Note that the Q x load is reacted by the facesheets as 
well as through the web, resulting in much lower stresses. At the same time, the Q y load is reacted strictly 
in the facesheet meaning that the x yz stresses are much higher. Figure 24 provides a detail of the stresses in 
the facesheet/flange combination, still in global coordinates. Figure 25 shows ply stresses in ply local 
coordinates, tn and T 23 . These stresses are used to calculate margins of safety, which are plotted in figure 
26 for the facesheet. Note that the margin of safety in the I-stiffener have not been plotted as they are very 
high. The minimum margin of safety is 1.74, indicating that the panel is over designed for the considered 
loading. 



Figure 22. — Location of the cut A-A in the l-stiffened panel. 
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Figure 23. — Shear stress distribution in the l-stiffened panel at 
cut A-A referred to panel global coordinates. 



Figure 24. — Detail of the shear stress distribution in the l-stiffened 
panel at cut A-A referred to panel global coordinates. 
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Figure 25. — Detail of the shear stress distribution in the l-stiffened 
panel at cut A-A referred to local ply coordinates. 
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Margin of Safety 

Figure 26. — Margin of safety distribution within the l-stiffened panel facesheet. 

Design Example 4 — Unstiffened Laminates in Pressure Shell Structure 

This example investigates the use of the simplified shear solution in the design of a complicated 
structure, in which the loads applied on each component are obtained from a finite element analysis. 
Subsequently, the structural sizing code HyperSizer (Collier Research Corp., 2007) is used to interface 
with the finite element model of the structure, extract the appropriate component loads, and then size the 
components based on obtaining all positive margins of safety. The simplified shear solution has been 
implemented within HyperSizer to provide an additional ply level margin of safety (see eq. (19)). 

The structural components considered are the 14 unstiffened laminate panels within a spacecraft 
pressure shell. The composite material considered is a woven graphite/epoxy laminate, whose ply 
material properties, which are needed by both HyperSizer and the simplified shear solution, are given in 
table 4. The ply thickness is 0.008 in., and HyperSizer is used to determine the number of plies (and thus 
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total thickness) needed for each composite component. Furthermore candidate laminates consisting of 
only 0° and 45° plies in a symmetric layup were considered. 


TABLE 4.— WOVEN GRAPHITE/EPOXY PLY PROPERTIES 


Eu 

Msi 

e 2 , 

Msi 

V12 

G n , 

Msi 

Ftu\, 

ksi 

Ftu 2 , 

ksi 

Fsu n , 

ksi 

Fsui 3, 
ksi 

Fsu 23 , 

ksi 

10.5 

10.5 

0.036 

0.84 

59 

59 

13 

6.9 

6.9 


Figure 27 shows the thickness of each component as obtained from FlyperSizer without consideration 
of the simplified shear solution. That is, no check has been made for ply failure vs. the through-thickness 
shear failure mode. It can be readily seen that most components have sized to a thic kn ess of between 0. 16 
in. and 0.416 in. (20 to 52 plies). Three components require a thic kn ess of approximately 0.9 in. 

The resulting simplified shear solution shear stress field for one particular component (as indicated by 
the arrow in fig. 27) are shown in figure 28. As can be observed, the x v - component is higher than the i xz 
component, which is due to the higher Q y load value for this component that was determined by the finite 
element solution. As indicated by table 4, the allowable shear stress (in ply coordinates) is significantly 
lower (6.9 ksi) than the maximum shear stresses magnitude indicated in figure 27 (51.14 ksi). Thus, 
failure of this component can be expected due to the high value of the induced shear stress. As shown in 
figure 29, when the simplified shear solution failure mode is considered, 1 1 of the 14 components of the 
structure are predicted to fail, as indicated by negative margins of safety (red color) in figure 29. 

Next FlyperSizer is permitted to size all laminate components of the structure vs. the simplified shear 
solution failure mode (in conjunction with all other failure modes). The resulting laminate thicknesses are 
shown in figure 30. This figure shows that much greater thicknesses are required for a closed design (i.e., 
no negative margins) when the simplified shear solution failure mode is considered (compare to fig. 27). 
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Figure 27. — Unstiffened laminate component thickness (in.), 
sized without simplified shear solution failure mode. 
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Figure 28. — Shear stress distribution within the component indicated by the arrow in figure 27. 



Figure 29. — Margin of safety of the previously designed composite laminate components. 
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Figure 30. — Unstiffened laminate component thickness (in.), 
sized including simplified shear solution failure mode. 


Conclusion 

A simplified theory based on the strength of materials approach has been presented for determination 
of the shear stress distribution in laminated plates and panels caused by applied shear force resultants. The 
method is straightforward and does not require the solution of elasticity equations in the presence of the 
corresponding boundary conditions. As such it provides an estimate of the shear stress distribution with 
no knowledge of the details of the structure using only a panel cross-section geometry, material 
properties, and the applied load. 

The predicted shear stress distribution for laminates are shown to be in excellent agreement with two 
alternative more complicated solutions, namely HOTFGM and HyperSizer Joints, for isotropic plies. For 
anisotropic plies, the agreement between the simplified shear solution, FlyperSizer Joints, and the exact 
(Williams, 1999) is satisfactory given the simplicity and efficiency of presented expressions involved in 
this method. 

Finally, in four design examples, the method has been applied to predict the biaxial shear distributions 
in an unstiffened laminate, a honeycomb sandwich panel, and I-stiffened panel. The shear stress 
distributions in the local ply coordinates are obtained by the relevant transformations, enabling 
determination of the margins of safety of these composite structures. 
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Appendix A — Plane Strain Effective Young’s Modulus Replacement 

The fully multi-axial orthotropic normal constitutive equations can be written as, 

(21) 
( 22 ) 
(23) 

In both beam theory and plate theory, the through-thickness stress, a., is assumed to be zero, allowing the 
eqs. (21) and (22) to be written as, 



-a r - 


v yx 


(24) 


v xy 


G x +- 


(25) 


In beam theory the y-dircction assumption is plane stress, ct v = 0, leaving for eq. (24), 


(26) 


In plate theory, the y-dircction assumption is plane strain, s v = 0, allowing a v to be determined from eq. 
(25) as, 


n v -t>’ , 1 

0 = — -^g t + — 


-» 


" v xy a x 


(27) 


Substituting eq. (27) into eq. (24) gives, 


-cJv - 


' yx 


~ v x y a x 


1 V XV V yx 

-> £ x = 


(28) 


Comparing eq. (28) to eq. (26), we see that the correct substitution in going from plane stress (beam 
theory) to plane strain (plate theory) is, 


1 V xyV yx 


(29) 


Also, noting that, 
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( 30 ) 



xy 


we have, 


E x 



( 31 ) 
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Appendix B — Shear Sign Convention 


The sign convention for positive external shear loading is shown in figure B.l. This is consistent with 

eq. (4), 


ct * = 


Pxz 

I 


which indicates negative (compressive) normal stress in the beam for a positive shear load (P) in the 
lower half of the beam (z < 0). The sign convention for the internal shear force (V) in the beam, which 
corresponds to the sign convention employed by Beer and Johnston (1981), is shown in figure B.2. As 
shown, the reaction shear force ( R ) at x = L is equal to P, while, due to equilibrium, the internal shear 
force in the beam, V= — P. Taking a horizontal cut of one of the portions of the beam, it is clear that the 
indicated sign convention for the directionality of the internal horizontal shear force (H) is consistent with 
the vertical shear force (V) sign convention. Considering an infinitesimal element of the beam along the 
face of the horizontal cut, it is clear that the standard sign convention for the internal shear stress is 
opposite of that employed for the internal shear forces. This is the motivation for the negative sign in 
eq. (11), which then enables a positive applied shear force P to result in a positive internal shear stress 
distribution. 


z 



Figure B.l . — Positive sign convention for external shear loading on a beam. 
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Figure B.2. — Positive sign convention for internal shear force and shear stress in a beam. 
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